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O^ ' We formulate the non-linear field theory for a fluctuating counter-ion distribution in the pres- 

^\ , ence of a fixed, arbitrary charge distribution. The Poisson-Boltzmann equation is obtained as the 

^\ • saddle-point, and the effects of fluctuations and correlations are included by a loop-wise expansion 

around this saddle point. We show that the Poisson equation is obeyed at each order in the loop 
^ , expansion and explicitly give the expansion of the Gibbs potential up to two loops. We then apply 

q_J ' our formalism to the case of an impenetrable, charged wall, and obtain the fluctuation corrections 

[L I to the electrostatic potential and counter-ion density to one-loop order without further approxi- 

mations. The relative importance of fluctuation corrections is controlled by a single parameter, 
which is proportional to the cube of the counter-ion valency and to the surface charge density. We 
also calculate effective interactions between charged particles, which reflect counter-ion correlation 
f 1 ^ ■ effects. 

o 

^ . PACS. 82.70.-y - Disperse Systems. 
^ ■ PACS. 61.20.Qg - Associated Liquids. 
PACS. 82.45.-Fz - Electrochemistry. 
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I. INTRODUCTION 



The behavior of charged, fluctuating systems is an old problem in chemistry and physics and is of importance for 
very diverse disciplines. In this article we will be concerned with the distribution of counter ions around charged 
iy~N objects, which is experimentally realized whenever an object with dissociable surface groups is brought in contact 

OO with water (or some other high-dielectric solvent). In such a situation, the counter ions will be attracted to the 
^^ i oppositely charged object, but also repelled from other counter ions. The characteristic feature of charged systems is 
^^ ■ that the Coulomb interaction is long-ranged, which gives these systems very special properties. Since, in general, one 
^r^ counter ion interacts with many different counter ions simultaneously, the mean-fleld approach is very successful and 
Qs ' can be used to describe experiments or simulations quantitatively. Exact solutions of the so-called Poisson Boltzmann 
"^ . equation, which determines the electrostatic potential distribution on a mean-fleld level, are available for planar |3^q] 
C^ ' and cylindrical U] geometries. A very readable introduction into the Poisson-Boltzmann (PB) approach is given in 
H Ref. [g. There are several factors which contribute to deviations from the PB equation, including additional, short- 

ranged interactions and solvent effects. In this article we will consider deviations due to fluctuations and correlations. 
These effects have been incorporated for the single charged wall by modified PB equations [|6|, by integral equation 
Q theories [[7|,y|, and by numerical methods |]9Hl2|. In general, one finds correlation effects to become important for 
^ ions of high valency, and packing effects dominate the ion distribution for very bulky ions or high ion concentration. 
L| For the case of counter ions or electrolyte solutions between two charged walls, there have been a number of studies 
• I— I , using integral equations ||l3|-[l5| and Monte-Carlo methods ||l^,|l3l ■ An exhaustive review has been given by Attard 
^% [Q . The main result of these studies is that correlations can lead to an effective attraction between similarly charged 

n3 : walls for high enough ion valency. There have been a few attempts of using field theoretic methods to go beyond 
■ ■ ' mean-field theory for charged systems. Here we only mention the works by Podgornik and Zeks [ p9[ and by Attard 
and coworkers pC| , who expanded around the mean- field solution and calculated within additional approximations 
fiuctuation corrections to the interaction between two charged, planar walls. 

In this article we consider a system of mobile, point-like counter ions in the presence of a fixed charge distribution, 
and formulate the field-theoretic framework for going beyond the Poisson-Boltzmann description. As has been realized 
by many others before, the Poisson-Boltzmann equation constitutes the saddle point of the exact field theory. We 
first reformulate this theory by a loop- wise expansion around the saddle point. The Poisson equation, which relates 
the electrostatic potential to the counter ion distribution, is obeyed at each level of this loop expansion. Clearly, 
the counter- ions are not distributed according to the Boltzmann weight, and correlations and fluctuations lead to 
pronounced deviations from the Poisson-Boltzmann predictions. The description is considerably facilitated by a 



Legendre transformation, after which we obtain the Gibbs potential as a function of a fixed electrostatic potential 
distribution. This transformation restricts the number of diagrams to so-called one-particle-irreducible diagrams. The 
electrostatic potential is given by the solution of the equation of state. 

We then apply this general formalism to the case of counter ions confined to a half space and next to a charged wall. 
We calculate without further approximations the one-loop correction to the electrostatic potential distribution and the 
counter-ion distribution. Due to correlations the counter ions are more densely packed close to the wall than predicted 
by the PB solution. This is in accord with previous Monte-Carlo simulations and with the intuitive expectation, since 
mean field usually overestimates repulsive interactions due to the neglect of correlations. The relative importance 
of fluctuation corrections to the electrostatic potential is measured by the single parameter q^a(.%, where a is the 
surface charge density of the planar surface, £b is the distance at which two elementary charges interact with thermal 
energy, and q denotes the counter-ion valency. Since the thickness of the counter ion distribution is, according to 
the mean-field prediction, given by /i ^ l/^qaiB), the average counter- ion concentration follows as c ^ f/^A* ^^^d 
our effective fluctuation parameter can be rewritten as q^cri% ~ {q'^iBC^'^)^'^ and thus measures the electrostatic 
interaction energy between two counter ions at their average separation (in units of the thermal energy); it is thus 
related to the plasma parameter. For vanishing values of q^crt"^ mean field theory becomes exact and fluctuation 
corrections are unimportant. On the other hand, since the counter ion valency comes in as a cubic power, it is 
clear that by going from monovalent to divalent ions, fluctuation effects become much more pronounced, in accord 
with experiments and simulations. Our results, which are obtained only to first order in the loop expansion, become 
unreliable for large values of q^a£^. A clear break down of our expansion is indicated by an unphysical negative 
density of counter ions, which occurs at q'^ai^ ~ 12. 

In Section II we formulate the general field-theoretic description for an ensemble of fluctuating counter ions in the 
presence of a fixed charge distribution. The main steps consist of i) formulating the initial partition function as a field 
theory using a Hubbard-Stratonovich transformation, ii) going to the grand-canonical ensemble, iii) expanding the 
action around the saddle-point, and iv) performing another Legendre transformation after which the Gibbs potential 
depends on the electrostatic potential. In Section III we show how the electrostatic potential follows from the equation 
of state, and in Section IV we give explicit results for the case of a single charged wall. In Section V we calculate the 
effective interaction between charged test particles, and Section VI is devoted to a brief discussion. 

II. NON-LINEAR FIELD THEORY FOR CHARGED SYSTEMS 

The partition function for N mobile counter ions of valency q, interacting with an arbitrary fixed charge distribution, 
(T(r), reads 



^^-M 



Jl / dr, cxpl-q^J2'"i^J-^k)+qdra{r)J2Hr~^j)}- (1) 



We assume the counter ions to be point-like particles, and only include the Coulomb interaction between ions, 
v{r) = £b/t, where £b is the Bjerrum length defined as £b = e^ /{A-KekBT). We also neglect any image charge effects 
in the present formulation of our theory. The charge distribution cr(r) in general also contains test particles, the 
presence of which allows to calculate the effective interaction between charges, as we will demonstrate in Section VI. 
Introducing the particle density operator 

N 

p(r)^^<5(r-r,), (2) 

we can rewrite the partition function as 

N 



ZN[h] = ^ 



n dr 



exp 



[-- I dvAr'[qp{v) - a{Y)]v{v - Y')[qp{v') - a{Y')]+ f drp{r)h{r)\ . (3) 



There is no need to subtract the diagonal terms involving self interactions at this point, because they will be canceled 
by a one-point rcnormalization at a later stage. Using the generating field /i, the average counter ion density can be 
calculated as 



(/5(r)> 



SZNJh] 
ZN[h]6h{r) 



(4) 



h=0 



Correlation functions can be calculated by taking multiple functional derivatives with respect to h. Noting that the 
operator-inverse of the Coulomb interaction is 



?;-i(r) 



47rf R '' 



the partition function is after a Hubbard-Stratonovich transformation given by 



I f V(J3 , 



dr 



S<'«' 



zcr(r)(?!)(r) 



^YQ-^'l>(^)+h{T^) 



N 



(5) 



where we introduced the rescaled Bjerrum length £b = A-Kq^ls- In the above expression, Zq denotes the partition 
function of the Coulomb operator, Zq ^ Vdet v. Using the definition Eq. (^, the average particle density turns out 
to be 



(p(r)> = N 



-l<t>{T) 



/dre- 



(6) 



and the normalization property of the density distribution, J drp{r) = N, is self-evident. The two-point cumulant 
correlation function (for r ^ r') follows as 



(p(r)p(r')) 



S^ZN[h] 



ZN[h]Sh{r)6h{r') 



/i=0 



N{N - 1) 



-t(p{r) — i(p{r') 



[Jdre 



-i4>{r)] 



(7) 



and higher-order correlation function can be calculated in a similar manner. The functional form of the particle 
distribution, Eq.(|6|), shows that icj) is the reduced electric potential. The partition function is brought into a more 
manageable form by going to the grand-canonical ensemble. 



Z[g,h] = y X^ZN[h] = [^ e^p\~in[(l),h] + ii /dr(/)( 
i^o J Zo I J 

where A is the fugacity of the counter-ions, and the Hamiltonian is defined as 

n[c^,h] ^ /dr [^(V^)^ + !^i^M£) _ Ae-^^«+"w' 
J [2£b q 



r)g(r) 



(8) 



(9) 



and we added a source term to the partition function in such a way that one can calculate directly the reduced 
electrostatic potential 



^(r) = (z0(r)) = 



5 In Z[g, h] 
£5g{r) 



(10) 



In the grand-canonical partition function, Eq.(|8|), we arbitrarily multiplied the action by a constant i which plays 
the role similar to the inverse Planck's constant in front of the action integral in Quantum Field Theory. This 
constant serves as an expansion parameter in our systematic treatment of fluctuation effects and counts the number 
of diagrammatic loops. The fugacity A is related to the total number of ions by 



A 



a In Z[g,h] 
IdX 



= (N) = X dr /e- 



(11) 



which establishes a useful relation between A and N. The particle density is obtained from the grand-canonical 
partition function Eq.(Q) using the definition Eq.(Q) and reads 



(p(r))=A(e-^W). 



Likewise, the two-point correlation function (for r ^ r') reads 

(p(r)p(r')>=A2(e-»'^W-.0(r')\^^ 



(12) 



which is simpler than the canonical form, Eq.(0); higher-order correlation function can be calculated similarly. In 
the following we will be interested in calculating the Gibbs potential r[tp,h] which depends on the reduced electric 
potential V'(r); it is related to hiZ[g,h] by a Legendre transform 



r[V', h] + In Z[g, h]=£ dr 'ip{r)g{r). 
Using the definition (llfl) we obtain from the definition of the Legendre transform (O) the inverse relation 



9ir) 



£Si^{r) 



(13) 



(14) 



We will perform the Legendre transformation using a loop expansion, using methods developed by Schwinger |2l| ] 
and closely follow the notation of Brezin, Le Guillou, and Zinn- Justin ||22(| . To that end, we expand the Hamiltonian 
around a fixed value of the potential, the so-called classical field, denoted as V'ci(r)- At this point, we do not specify 
what this classical field is. We write 



(r) ^ -i->pci{r) + xir) 



(15) 



and thus obtain the formal expansion 



n[^] ~ij Av<i>(v)e{v) = H[-iM - jdrgir) [V^ei(r) + *x(r)] + ^ }-j'H^J){{v.^}- ^,{)x{ri) ■ ■ ■ x(r,)dri • • • dr,. 
The vertex functions Ti'-') are defined by 



HW({r,};V'ci)^ 



5^1-L\ 



<50(ri)...(50(r,) 



(16) 



6— — ^^c 



and of course depend on the function V'ci(r) around which one expands. The function lnZ[gi] can now be evaluated 
by standard saddle-point methods, treating l^^ as the expansion parameter. The saddle point i/'sp(r) is defined by 



5n\ 



(r) 



ig{T) 



= 



(17) 



0=-j0sp[e] 



and depends on the source g. Using the definition of the two-point propagator G(r, r'; V'sp); 

j dr G(r, r'; Vsp)W<'^ (r, r"; ^gp) = ^(r' - r"), 
we can diagrammatically expand the expression for hiZ[g, h] and obtain up to two loops 

lnZ[g,h]=l /"dr(^(V^sp)'-^^^^^^^^^^ + Ae-^^^W+''('-)-fe(r)i/;sp(r)l-ilndet?^(2)[^gp] 

i |-i y dn • • • dr4 n^^Kiuh Vsp)G(ri, r2; Vsp)G(r3, u; V'sp) 
+ y dri • • -drgdr; • • • dr^ H(3)({r3}; Vsp)H(3)({r^}; V^Sp) 

-G(ri,r2;Vsp)G(r;,r'2;Vsp)G(r3,r^;^Sp)-)-— G(ri,r;;^Sp)G(r2,r^;7/>sp)G(r3,r^;7/>sp) 



-I- 



(18) 



The one- loop diagram and the two- loop diagrams are schematically represented in Fig.l. In order to perform the 
Legendre transformation, we need the loop expansions for V'(r) and g{r), which are obtained from the definition (M) 
and the saddle-point equation (p7|), respectively. These expressions depend on the saddle-point potential ipsp- We 



then have to expand all functions and vertices around the saddle point and reexpress the potential dependence as a 
function of the expectation value f/;. By inserting these expression into the Legendre transform, we obtain the loop 
expansion of the Gibbs potential, 

V[^,h] = -I fdrl^iV^pf - ^ilMll + Ae-^W+'^MJ + ^ IndetT^^^)^^] 

-i|-iydri...dr4W(^)({r4};^)G(ri,r2;^)G(r3,r4;V) 

+ _L y dn • • • dradr'i • • • dr^ n^'Kivs}; ^/jW\{r',}; i>) 
G(ri, ri; V')G(r2, r'2; ^)G(r3, r^; V)} (19) 

The effect of the Legendre transformation is to cancel all one-particle-reducible diagrams, as expected. For the present 
case, the two- loop diagram to the right in Fig.l is removed from the expansion. 

III. THE EQUATION OF STATE 

The equation of state is defined as 

and completely determines the electrostatic potential. In the following, we will restrict ourselves to the one-loop order. 
Using the explicit form of the two-point vertex function, which follows from the definitions Eq. (O) and Eq. (Uq) , 

7^(2)[i.,r';V]= ('-^ + Ae-'^(""))(5(r-r'), (21) 



the equation of state can be explicitly written as 

- Ae-^W - ^ - ^^M + Ac-VHr)G(r, r) = 0. (22) 

q Ig 2£ 

Note that from now one, we suppress the dependence of the two-point correlation function G(r, r') on the electrostatic 
potential distribution ip. The correlation function is determined by the equation 

-^-f Ae-'^('"))G(r,r') = <5(r~r'). (23) 

^B J 

Finally, the number of mobile ions can be calculated according to Eq. ( pi] ) and is given by 

{N)^\ I dr e-'^f'') _ A /" dr e-'''('")G(r, r). (24) 

The equation of state can be solved by a systematic expansion of all quantities which are to be determined in inverse 
powers of the loop parameter L At the one-loop order, only the leading term of G contributes and therefore an 
expansion of G itself is unnecessary. The electric potential and the fugacity have to be expanded to the next-leading 
order, 

V'(r)=Vo(r)-friVi(r), (25) 

A = Ao+riAi. (26) 

The equation of state splits into two separate equations, the zero- loop equation 

- Aoe-'^»(r' - ^ - ^!|^ = 0, (27) 

q iB 



which is the ordinary Poisson Boltzmann equation, and the next-leading correction, 

_ x,e-^o{r) ^ Aoe-^°Wv^i(r) - Z!^ + ^e-^°«G(r,r) = 0. (28) 

The correlation function is determined by 

72 



(-Z- + Xoe-^^^'A G(r, r') ^ S{r - r'). 



(29) 



The fugacity is at the zero-loop level determined by 



^° /dre-^o(r) (30) 



and at the next-leading level by 



^ _^ /dre-'>°WV;i(r) , Ap /dr e-^"WG(r,r) 

^^"^" /dre-^o(r) +2 /dre-'AoW " ^"^^^ 

Combining the equations for -01 and G, Eqs.(p8|) and (E9), we obtain for the correction to the electrostatic potential 

2Ar 



01 (r) = -:^ / dr'e-^''(''''G(r',r) 



G(r',r') 

-^0 



(32) 



which contains an implicit dependence on ipi on the right-hand side through the dependence of Ai on ipi, see Eq.(pl|). 
Algebraically solving Eqs.(Bl]) and ( p2| ) we obtain 

dr e-'>o(r) , . . ^ ~(Ao/2)/drdr' e-'>»W^'>°C-')G(r-,r) [G(r,r) - Z/F] 

^'^ ^ l-(Ao/F)/drdr'e-^o(r)-^o(r')G(r',r) ' ^ ' 

where we used the short-hand notations 

Y= /"dre~'^°('") (34) 

and 

Z= /"dre-'^''WG(r,r). (35) 

The one- loop correction to the electrostatic potential is therefore determined by Eq.(p3) in conjunction with the 
equations determining G and Ai, Eqs.(p9[) and (O). 



IV. SINGLE CHARGED IMPENETRABLE WALL 

In the following we consider the case of a charged wall which is impenetrable for the counterions. The fixed charge 
distribution is given by a{r) = a 5{z) where a is the surface charge density at the charged wall. The solution for the 
one-dimensional Poisson-Boltzmann equation, Eq. ( |27| ) , is 

0o(^)=21n(l + z^), (36) 



where /i = \^2/tB^o = l/-\/27rq2^sAo is the Gouy-Chapman length. Enforcement of the normalization condition, 
Eq. (pO[) , leads to a relation between the surface charge density a and the Gouy-Chapman length, a = qXaH, which 
can be further transformed to Ag = 2T:(.Ba^ and l//x = 2T:iBqu, the standard results. 



A. Calculation of the correlation function 



Due to lateral invariance, the Green's function G(r, r') can be Fourier transformed in the directions perpendicular 
to the surface; the equation for the Green's function, Eq.(p9[), becomes 

dz2 



P 



Xoe 



CBAo 



^i'oiz) 



nG{z,z',p)^£BS{z~z'). 



(37) 



where p are the Fourier variables in the lateral directions, and p = |p|. The potential ipo{z) is given by Eq. (p6[) for 
z > 0. Since the ions cannot penetrate into the wall, the counter-ion density is zero for z < 0, or, equivalently, we 
put Tpoiz) — oo for z < 0. The solution can be obtained in a straightforward manner and reads 



G{z,z',p) = 



P 



2p5 \^ z' + ^ 



P 



P 



M 



p^ I e 



P 



.p(z^z') _(_ I Z + 11 ^ 



p" 



^-p(z+z') 



(38) 



for z' > z > 0. The solution for z > z' > follows from Eq. 
limit p = we obtain 



1 + 2p2^2 _^ 2p^ 
by an interchange of the arguments z and z'. In the 



G(z,z',p = 0) = 



lB^i{'2 + {l + z/^if) 



(39) 



3(l + z//i)(l + z'/Ai) 

for z' > z > (and with the coordinates interchanged for z > z' > 0). The self interaction is defined as the equal-point 
Green's function and follows from the partially Fourier-transformed Green's function as 

dPp 



Gsc\i{z) 



{M 



:G{z,z,p). 



The integral to be solved is, subtracting the (infinite) Coulomb self energy Gsc\i{z ~ oo), 
^ ,, Ib r dp \ 1 , l + 2p-\z + ^i)-^+p-\z + ^i)-^ 







27r 



-2pz 



(40) 



(41) 



(z + fi)'^p'^ 1 + 2p2/i2 -|- 2pfi 

The integral can be performed by some intermediate transformations, and the result can be written in a scaling form 
as 



tB 



Gsciiiz) = -^g{z/fi) 



where the scaling function g(x) is explicitly given by 

1 



9ix) 



ie^^-'i}^Ei[{l - i)x]{l + ixf - ie(i+')^Si[(l + i)x\{l - ixf - 4x 



(42) 



(43) 



2(l + x)2 

and plotted in Fig. 2. The function Ei[x] is the exponential-integral function EJ. The prefactor of the scaling function 
in Eq.(|4|), 

27rq3£|cr, (44) 



47r/^ 



M 



is a measure of the importance of correlation effects and thus the departure from mean-field behavior. We note that 
it is this single parameter combination which measures the relative importance of fluctuation corrections. Clearly, 
the cubic dependence on the counter-ion valency shows that fluctuation effects will become much more dominant for 
multivalent ions. As can be seen from Fig. 2, the self energy is negative, which reflects the fact that the screening due 
to counter ions lowers the local free energy of charged particles. The scaling function g{z/ n) shows a minimum at 
a finite distance from the wall, z/^ ~ 0.2. This maximum in screening is closely connected with a maximum in the 
counter- ion density at about this distance from the wall, as we will show in Section IV. C. The asymptotic behavior 
for small separations from the wall is 



9{3c) 



4 ^ ' 



and 



/W 



3 

2x 



(45) 
(46) 



for large separations. The functions Eqs.(^5|) and (E^) are shown in Figs. 2a and 2b as broken lines. 



B. Calculation of the one-loop correction to the electrostatic potential 

The solution of the equation for ipi is considerably simplified by the fact that, for the specific mean-field potential 
given in Eq.(|3^), we find 

dz G{z, z',p ^ 0)e-'^°(") = ^~bmV2 (47) 

independent of z'. It immediately follows from Eq.(p3) that the integral /dze~'^°^^Vi(^) is zero and, as a consequence, 
the equation for -01, Eq.(p3), is independent of any constant appearing in the correlation function G. This holds 
in particular for the Coulombic self energy and therefore justifies that we neglected this (infinite) constant in the 
calculation of the self interaction. This is an example of a renormalization of a one-point correlation function, which 
otherwise is plagued by the (infinite) Coulombic self-energy. The normalization of the bare potential V'o i Eq. ( |36| ) is 

dz e~^°'^'^ = fi. (48) 



(49) 



The renormalization constant appearing in Eq. (|32D thus is 

2Ai ^ Jdze-^o(^^GMz) ^ JBCi 
Ao /dze-'^ol^) 47r^ 



where the numerical constant is ci ~ L dxg{x){l + x) '^ ~ 0.6208593. Our final result for the potential is, according 
to Eq.dl), 

'^i(^) = ^[/(^/M)-y], (50) 

where the scaling function f{x) is defined by 

^^"^ ^ " X '" m+x)[i+x'Y '^^ ) - 1 ^^ 3(i+x)(i+.o^ g^" ) (^^) 

and graphically presented in Fig. 3. The asymptotic behavior for small separations is 

/(.)-3(l-^ (52) 

and for large separations we obtain 

f{x) c. i^. (53) 

The asymptotic formulas are shown in Fig. 3 as broken lines. 

C. Calculation of the counter-ion density 



The counter ion density p(r) can be calculated from the Gibbs potential Eq.(|19[) by taking a derivative with respect 
to the generating field h. The result is, including the one-loop correction, given by 

p(r)=Ae-^W-Ae-'AMG(r,r) (54) 



and leads via integration to the total ion number TV, Eq.(g4|), as it should. Subtracting the equation of state, Eq.(P2[), 
we obtain the Poisson equation, 

p(r) = -^-X!Vl£). (55) 



In fact, the equation Eq.(|5^) is exact to all orders in the loop expansion, since the equation of state contains exactly 
the same terms as does the particle density and thus leads to an exact cancelation of all higher-loop terms. This is 
so because in the Hamiltonian Eq.(|^), the generating field h enters in the same way as the fluctuating potential </>, 
except for the linear and quadratic terms. Expanding the counter ion density in powers of the loop parameter, 

p(r)-po(r)+rVi(r), (56) 

and using the loop-wise expansions of the fugacity A and the electrostatic potential ■(/'(r) introduced in Section III, 
the zero-loop (or mean-field) result for the counter ion density is 

Po(r) = Aoc-'^"('-) (57) 

and the one-loop correction reads 



pi(r) = Aoe-^°M ("^ - Vi(r) - \g{v, v)\ 



(58) 



Using the explicit solution for the single charged wall, Eq.(pfl), and the result for the mean-field fugacity, Aq = IjlBli^ ■, 
the zero-loop density can be written as 



Po(^) = ( y^) ,,2 I ,, (59) 



and the one-loop correction reads 



where the scaling function h[x) can be written in terms of the previously defined scaling functions and reads 

In Fig. 4 we plot the scaling function h{x). The density correction vanishes at the charged wall, and for small 
separations the asymptotic behavior follows from Eqs. (^5|) and (p^ ) as 

TT 

h{x) ~ — — a;lna;. (62) 

8 

The density corrections is positive for x < 1, i.e. for distances from the wall smaller than the Gouy-Chapman length 
/i fluctuations and correlations lead to an enhanced density, quite in accord with expectations: Each counter ion is 
surrounded by a correlation hole, which is neglected on a mean-field level. Mean field theory therefore overestimates 
the repulsion between counter ions, and therefore underestimates the density close to the charged wall. This is in 
qualitative agreement with Monte-Carlo simulations [Q. We note that since the slope of the density correction is 
infinite close to the wall, see Eq.(|62|), our results predict always an initial increase of the counter-ion density as one 
moves away from the substrate. The total integral over the density correction vanishes; the increase in density close 
to the wall is therefore compensated by a decrease in density further away from the wall, as shown in Fig. 4b. For 
large separations from the wall, the asymptotic behavior follows from Eqs. ( p5| ) and (|2|) as 

h{x) c -i^. (63) 

x'' 

In Fig. 5 we plot the rescaled density p{z) = (i?B/i^/2)p(z) as a function of the rescaled distance from the wall z/fi, 
for four different values of the parameter ^B/47r/i. The broken line, for i?B/47r/i = 0, denotes the PB result, the other 
three solid lines are for £B/47r/i = 1, 5, and 10, with the difference to the broken line increasing as the parameter 
IbI^t^I^ increases. We see that already for rather small values of £3/477/1, the one-loop result for the counter ion 
distribution is quite different from the PB result. The density depression far away from the wall is proportional to 
the parameter £3/477/1. It is clear that for very large values of this parameter, the sum of the zero-loop and the 
one-loop densities will become negative. This in fact happens at £3/477/1 w 12, which clearly sets an upper limit to 
the applicability of the present one-loop calculation. It is likely that the one-loop calculation becomes quantitatively 
inaccurate even for smaller values of this parameter, but the precise value of £3/477/4 when this happens can only be 
judged by a numerical calculation (such as Monte Carlo) or by a higher-loop calculation. 



V. INTERACTION BETWEEN CHARGED PARTICLES 

In this section we calculate the effective interaction between charged test particles. As mentioned before, test 
particles can be incorporated into the model by adding them to the charge distribution in Eq.([l|) according to 

M 

a{r)^a{r) + J2QjS{r-Rj), (64) 

where we have considered a collection oi j = 1, . . . ,M test particles with charges Qj and positions R^. Inserting this 
shifted charge distribution into the grand-canonical partition function Eq. (^), the reduced test-particle free energy 
turns out to be 

/({Rm}) = - hi /c^^S. Q.^(R.)/A . (65) 



Using the separation of the fluctuating field (j) into the expectation value ip and the fluctuations around the expectation 
value, Eq.([lq), the free energy can be reexpressed as 



M -, M 



/{{Rm}) = J2 QM'^j)/'! + IJ2 Q,QkGiR,,Rk)/q^. (66) 

Using the loop-expansion of the electrostatic potential, tp, Eq. (p5|) , and the definition of the self interaction, GgcU, 
Eq. (pO[) , the free energy becomes on the one- loop level 

/({Rm}) - E j— [MZj)+MZj)] + ||G.eif(Z,)| +i^ ^G(R„R,) (67) 

and thus separates into a single-particle part and a two-particle part. The single-particle contribution contains the 



ordinary Gouy-Chapman potential, ipo{z), given by Eq.(36), a contribution due to fluctuation-induced changes of the 



electrostatic potential, ^i{z), given by Eq. (pO[) , and a contribution due to correlations between counterions and a test 



particle, Q^Gscii{z), given by Eq. ([42[) . Since the correlation contribution is attractive and proportional to the square 
of the test-particle charge, it will be dominant for large test-particle charge and lead to attraction even in the case of 
a similarly charged wall. This somewhat surprising behavior has in fact been seen in Monte-Carlo simulations [0|. 
We note that on the two- loop level, in addition a three-point interaction appears. 

Finally we will present results for the interaction between two particles in the neighborhood of the charged wall and 
thus under the influence of the loosely bound counter-ion cloud. To make the results somewhat more transparent, we 
will confine ourselves to particles which have the same vertical distance from the wall. The effective interaction then 
follows by Fourier transformation and reads 

G{z,r)= f^G{z,z,p)Jo{pr) (68) 



where the Green's function G(z, z,p) is determined by Eq.(|3^) and Jo{x) denotes the Bessel function of the first kind 
[P3| . For particles which are very close to each other, r < ii, and/or very far apart from the charged surface, z > r, 
the interaction is basically unscreened and given by 

G{z,r)~£Bq^/r. (69) 

For particles that are far apart from each other, r > /i, but close to the surface, z < ^, the effective interaction is 
given by 

G(.,.).^^fl-4i + 64 + ...). (70) 

The main feature is that the effective interaction is screened and thus reduced compared to the bare Coulomb 
interaction, but the screening is much weaker than for example in a salt solution and the interaction decays as 
an inverse cube of the distance. This behavior has previously been predicted by liquid-state calculations and is 
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interpreted as the effective dipole-dipole interaction between the test-charges and their associated counter-ion clouds 
p4|-p6[ . The effective interaction between two induced dipoles is unscreened because the lower half-space is assumed 
to be impenetrable to ions and thus allows the unattenuated passage of electric field lines. From Eq.([70|) it follows 
that the screening is maximal at a finite distance away from the wall. This is connected with the maximum in the 
counter-ion distribution, see Fig. 4. For particles that are far apart from each other, r > fj,, and also relatively far 
apart from the surface, fj, < z < r, the effective interaction is given by 

G(z,r)~— 3 (71) 

and thus increases gradually as one move away from the surface. At a separation from the surface which equals the 
intcrparticle distance, z ~ r, the effective interaction has the same magnitude as the bare Coulomb interaction, and 
we find the effective interaction to cross over smoothly to the bare interaction as given by Eq. ( |69|) . 

VI. DISCUSSION 

In this article we have formulated the non-linear field theory for mobile counter ions under the influence of fixed 
charge distributions. We showed explicitly that the Poisson-Boltzmann equation corresponds to the saddle-point 
equation, and how correlations and fluctuations can be accounted for by a loop-wise expansion of the action around 
this saddle point. We find the Poisson equation to be satisfied at each order in this expansion. Clearly, the Boltz- 
mann equation is not satisfied when going beyond the saddle point, and this is the reason for deviations from the 
Poisson-Boltzmann equation. Particularly useful is the introduction of the Gibbs potential r['0], which only contains 
one-particle irreducible diagrams (and thus reduces the number of diagrams to be considered) and which allows to 
directly calculate the electrostatic potential distribution via the equation of state. We applied our formalism to the 
case of a charged wall which is impenetrable to counter ions and obtained, to one-loop order and without further 
approximations, the electrostatic potential, the counter-ion distribution, and the effective interactions between test 
particles. We find that the counter- ion density is increased close to the wall as compared to the mean-field (Poisson- 
Boltzmann, PB) solution, however, right at the wall the density is unchanged as compared to PB. The increase close 
to the wall is due to correlations between counter ions, which are not captured in the PB approach and which reduce 
the repulsion between counter ions. The surface potential is increased at the surface. The relative strength of the 
one- loop correction is proportional to the single parameter Ib/'^tt^j, = 2Trq^£^a and thus is proportional to the cube 
of the counter ion valency. Our results become unphysical at ^B/47r/i « 12. The value of £b/4:Ttii up to which our 
results are accurate could only be inferred from quantitative comparisons with numerical simulations, which wc plan 
to do in the future. Our formalism can in principle applied to all situations where the PB equation can be solved in 
closed-form. It would be particularly useful to obtain the one-loop correction for a charged wall in a salt solution, 
because this would represent the bridge between Debye-Hiickel and Poisson-Boltzmann theory. 
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FIG. 1. a) One- loop and b) two-loop diagrams that enter the calculation of the logarithm of the partition function, Eq .(h8h . 
The two- loop diagram to the right is one-particle reducible and thus cancels out when going to the Gibbs potential, Eq.(h^ 

FIG. 2. Plot of the rescaled self-energy 3 as a function of the rescaled distance from the wall, z/fj,, as determined by Eqr(H2|). 
The asymptotic behavior for small and large wall distance, Eqs.(|45[) and mw, are drawn as broken lines in a) and b), respectively. 

FIG. 3. Plot of the rescaled correction to the electrostatic potential /, as defined by Eq.dM). The asymptotic results for 
small and large separation from the wall, Eqs.(|52[) and (p3), are denoted by broken lines. 

FIG. 4. Plot of the rescaled one-loop correction to the counter ion density h, as defined by the Eqs.(|60[) and (pll). One notes 
that the density change at the substrate surface vanishes identically. 

FIG. 5. Plot of the one-loop prediction for the rescaled density p{z) = [istJ-'^ /'2)p{z) as a function of the rescaled distance 
from the wall z/fi, for four different values of the parameter £B/4:TTfj.. The broken line denotes the PB result; the other three 
solid lines are for £b/4:-kjj. — 1, 5, and 10, with the distance to the PB curve increasing as £b/4:Ttp increases. 
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Figure 1, Netz and Orland 
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